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A  simple  method  to  invert  the  absolute  geostrophic  velocity  on  potential-density  sur¬ 
faces  is  developed  based  on  McDougall’s  (1988)  conservation  statements  on  poten¬ 
tial-density  surfaces  and  Chu’s  (1995)  and  McDougall’s  (1995)  P-vector  concept.  This 
method  has  a  capability  to  handle  a  system  with  diapycnal  velocities  and  non-coinci¬ 
dence  of  potential-density  and  neutral  surfaces.  The  validity  of  the  inversion  is  the 
existence  of  the  P-vector  spirals  instead  of  the  existence  of  the  velocity  spirals  (usu¬ 
ally  called  the  /1-spirals).  We  use  a  climatological  T,  S  dataset  for  the  Atlantic  Ocean 
from  the  National  Oceanographic  Data  Center  (NODC)  to  show  the  benefit  of  using 
P-vector  spirals  in  inverting  the  velocity  field. 
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1.  Introduction 

Our  understanding  of  the  mid-latitude  large-scale 
ocean  circulation  has  been  benefitted  by  a  set  of  papers 
by  Stommel  and  collaborators  (Stommel  and  Schott,  1977; 
Schott  and  Stommel,  1978;  Behringer  and  Stommel, 
1980),  Wunsch  and  collaborators  (Wunsch,  1978;  Wunsch 
and  Grant,  1982),  Killworth  (1986),  and  Chu  and  collabo¬ 
rators  (Chu,  1995;  Chu  et  al.,  1998a,  b).  Their  work  makes 
it  possible  to  obtain  ocean  general  circulations  from  ob¬ 
servations  of  temperature  ( T)  and  salinity  ( S ).  The  physi¬ 
cal  base  for  calculating  geostrophic  velocity  from 
hydrographic  data  is  the  thermal  wind  relation 
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where  ( x ,  y,  z )  are  the  three  axes  with  corresponding  unit 
vectors  (i,j,  k)  of  the  Cartesian  coordinates  pointing  east¬ 
ward,  northward,  and  upward,  respectively;  V  =  id/dx  + 
jd/dy ,  is  the  two  dimensional  gradient  operator;  V  = 
(n,  v),  VQ  -  (n0,  v0)  are  the  geostrophic  velocities  at  any 
depth  z  and  at  a  reference  depth  z0,  p  is  the  pressure,  g  is 
the  gravitational  acceleration,  and  /  is  the  Coriolis  pa¬ 
rameter,  which  is  a  function  of  latitude.  As  mentioned  by 
Olbers  et  al.  (1985),  the  quantities  T,  S  are  relatively  easy 
to  measure,  and  in  contrast  to  velocity  observations,  the 
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climatological  signal  in  the  T,  S  fields  is  less  contami¬ 
nated  by  energetic  smaller-scale  motions  induced  by  ed¬ 
dies  and  waves.  Usually,  the  hydrographic  data  only  de¬ 
termine  the  relative  geostrophic  currents  with  the 
Boussinesq  approximation. 
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The  reference  velocity  V0  still  needs  to  be  determined. 

Davis  (1978)  pointed  out  that  the  /1-spiral  method 
(Stommel  and  Schott,  1977)  and  the  box  method  (Wunsch, 
1978),  no  matter  how  different  in  appearance,  are  based 
on  the  same  order  of  dynamical  sophistication  and  differ 
from  implicit  assumptions  about  the  scales  of  oceanic 
variability  and  different  definitions  of  the  smooth  field 
to  which  the  dynamical  model  pertains.  The  physical  prin¬ 
ciple  for  both  methods  are  the  existence  of  a  conserva¬ 
tive  tracer  which  allows  determination  of  a  family  of 
material  (usually  the  potential  vorticity)  surfaces  z  = 
h{x ,  y)  such  that 


Nf'f'dz{V'+V0]  =  0.  (3) 


For  steady  flow  in  an  ideal  and  incompressible  ocean, 
density  (p),  and  potential  vorticity  ( q )  are  materially  con¬ 
served. 

Shortly  after  the  /1-spiral  method  was  proposed, 
Schott  and  Stommel  (1978)  reported  that  in  some  in- 
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stances  the  absolute  velocity  computed  at  a  constant  ref¬ 
erence  level  depends  on  the  depth  range  of  data  used  in 
the  computation.  Behringer  (1979)  thought  that  this  am¬ 
biguity  was  caused  by  large  error  in  calculating  0  in  the 
Z-level.  Here,  0  is  the  angle  between  the  horizontal  abso¬ 
lute  velocity  vector  and  the  eastward  direction.  He  pro¬ 
posed  a  two-step  procedure  to  reduce  this  ambiguity: 
(a)  determination  of  0  on  isopycnal  surfaces,  and  (b)  de¬ 
termination  of  the  reference  velocities.  Taking  the  varia¬ 
tion  of  (3)  with  respect  to  h  to  obtain  (Davis,  1978) 

V0'Vdzh-jV0  =  -dz(V'’Vh)jv'  (4) 

which  differs  from  the  original  jS-spiral  method  (Stommel 
and  Schott,  1977)  by  the  term  S7h’[d'V\.  Coats  (1981) 
used  (4)  to  calculate  the  absolute  velocity  in  the  north¬ 
eastern  Pacific  Ocean. 

Keeping  the  dynamical  system  as  geostrophic  and 
hydrostatic  balanced  and  adding  diapycnic  and  isopycnic 
turbulent  mixing  terms  in  the  tracer  equations  (i.e.,  heat, 
salt,  and  potential  vorticity  equations),  Olbers  etal.  (1985) 
applied  the  /1-spiral  technique  to  the  climatological 
hydrographic  atlas  of  the  North  Atlantic  Ocean  and  ob¬ 
tained  the  velocity  field  resemblance  of  the  classical  view 
of  the  North  Atlantic  circulation. 

The  intersections  of  isopycnals  with  surfaces  of  con¬ 
stant  q  must  parallel  to  streamlines,  and  thus  the  direc¬ 
tion  of  velocity  is  one  of  the  two  opposite  directions  (180° 
apart)  at  each  level.  Behringer  (1979)  appears  to  be  the 
first  to  employ  potential  density  and  potential  vorticity 
to  determine  flow  direction  but  did  not  combine  results 
from  different  levels  to  establish  velocity  (Davis,  1978). 
When  considering  the  absolute  horizontal  velocity  at  two 
levels,  the  magnitudes  can  be  determined  by  the  thermal 
wind  relation.  Apparently,  the  method  fails  if  isopycnals 
and  surfaces  of  constant  potential  vorticity  coincide  so 
that  the  velocity  direction  cannot  be  determined  (neces¬ 
sary  condition  1).  Another  case  of  failure  arrises  when 
the  absolute  horizontal  velocity  vector  does  not  turn  with 
depth  (necessary  condition  2)  and  so  incompatibility  with 
the  relative  velocity  profiles  can  occur  (Olbers  el  al., 
1985).  The  two  necessary  conditions  for  the  validity  of 
the  /5-spiral  method  can  be  easily  identified  by  a  unit  vec¬ 
tor 


of  necessary  condition  2.  Inversion  of  the  absolute  ve¬ 
locity  (V  =  yP)  is  divided  into  two  steps:  (a)  determina¬ 
tion  of  the  vector  P,  and  (b)  determination  of  y. 

McDougall  (1988)  developed  a  consistent  formula¬ 
tion  for  the  /5-spiral  in  an  isopycnal  framework.  Later  on, 
he  obtained  the  absolute  horizonal  velocity  on  the  neu¬ 
tral  surface  from  vertical  differentiation  of  the  potential 
vorticity  equation  with  the  thermal  wind  relation 
(McDougall,  1995).  Following  McDougall’s  (1988,  1995) 
framework,  we  develop  a  simple  method  to  invert  the 
absolute  velocity  on  isopycnal  surfaces  from  hydrographic 
data. 

The  outline  of  this  paper  is  as  follows.  We  introduce 
McDougall’s  (1988)  conservation  statements  on  poten¬ 
tial-density  surfaces  in  Section  2,  and  present  concepts 
of  the  P-vector  on  the  potential-density  surface,  two  nec¬ 
essary  conditions  for  the  inversion,  and  the  P-vector  spi¬ 
ral  in  Sections  3-5.  In  Section  6,  we  establish 
climatological  annual  mean  isopycnal  depth,  potential 
vorticity,  and  P-spiral  parameter  data  sets  for  the  Atlan¬ 
tic  Ocean  from  a  climatological  T,  S  data  produced  by  the 
National  Oceanographic  Data  Center  (NODC)  (Levitus 
and  Boyer,  1994;  Levitus  et  al.,  1994).  In  Section  7,  we 
investigate  the  inverted  latitudinal  volume  transport.  In 
Section  8  we  present  our  conclusions. 

2.  Conservation  Statements  on  Potential-Density 

Surfaces 

Consider  a  series  of  potential  density  ( ag )  surfaces, 
each  marked  by  a  constant  value  of  a,  and  with  the  depth 

z(a)  =  R( x,  y,  a). 

The  vertical  distance  between  two  closely-spaced  a  sur¬ 
faces  with  increment  of  Acr  is  given  by 


h(a)  = 


dz(a) 

da 


A  a. 


(5) 


The  conservation  on  the  potential-density  surface  is 
given  by  McDougall  (1988) 
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which  is  determined  by  the  density  field  only  (Chu,  1995; 
Chu  el  al.,  1998a,  b).  Existence  of  this  unit  vector  (P) 
implies  the  satisfaction  of  necessary  condition  1.  Exist¬ 
ence  of  the  vertical  turning  of  P  indicates  the  satisfaction 


where  Via)  and  V a  are  the  lateral  velocity  and  gradient 
operator  (two-dimensional)  along  the  potential  density 
surface;  Q(cT>  is  the  potential  vorticity;  qia)  is  a  conserva¬ 
tive  quantity  representing  the  potential  vorticity  (we  may 
call  it  the  pseudo  potential  vorticity),  and  w<0',  is  the 
diapycnal  velocity  through  a  (7-surface.  The  diapycnal 
velocity  w(<TI  has  contributions  not  only  from  the  vertical 
diffusivity,  thermobaricity  and  cabbeling,  but  also  from 
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lateral  mixing  along  the  neutral  tangent  plane  (McDougall, 
1988). 

McDougall  (1988)  further  derived  the  following  ther¬ 
mal  wind  relation  for  the  potential  density  surfaces. 


relation  (7).  Vertically  varying  diapycnal  velocity  (9w(<TV 
dz  =£  0)  causes  the  crossing  </,rf)-isoline  motion. 

Use  of  (6)  leads  to 
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where  Vis  the  buoyancy  frequency;  0is  the  potential  tem¬ 
perature,  and  n  is  the  factor  indicating  the  difference  be¬ 
tween  the  potential-density  and  neutral  surfaces.  The  neu¬ 
tral  surface  is  defined  as  a  surface  along  which  fluid  par¬ 
cels  can  be  moved  a  small  distance  without  experiencing 
a  gravitational  restoring  force.  If  fi  -  1,  the  neutral  sur¬ 
face  coincides  with  the  potential-density  surface 
(McDougall,  1988).  In  this  paper,  we  call  the  case  of  /i  = 
1  as  the  no  /i-effect  case  (no  difference  between  the  neu¬ 
tral  surface  coincides  with  the  potential-density  surface) 
and  the  case  of  fx  ^  1  as  the  /i-effect  case  (difference  be¬ 
tween  the  neutral  surface  coincides  with  the  potential- 
density  surface). 

3.  P- Vector  on  Potential-Density  Surface 

Following  Chu  (1995),  we  may  define  two  orthogonal 
vectors  on  the  potential  density  surface 
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where  Px  and  P  are  the  components  of  the  vector  P. 

Without  vertical  change  of  the  diapycnal  velocity 
equation  (9)  becomes, 


V(a)=yP  (11) 

and  any  water  particle  moves  along  r/a)-isolines  on  the 
potential  density,  i.e.,  any  gf<7'-isoline  is  a  trajectory  of 
water  particles  (Fig.  1).  We  call  this  case  (dw{a)/dz  =£  0)  as 
the  no  diapycnal  exchange  case. 

4.  Necessary  Conditions  for  Inversion 

Applying  the  thermal  wind  relation  (7)  to  any  two 
different  potential  density  levels  Ok  and  am  (the  corre¬ 
sponding  depths  zk  and  zm),  shown  as  in  Fig.  2(a),  a  set  of 
algebraic  equations  for  determining  the  parameter  y  is 
obtained 


which  are  the  tangential  and  normal  unit  vectors  of  the 
g^-isolines  (Fig.  1).  Using  (8),  the  velocity  vector  U^can 
be  decomposed  as 


f)p«)  _  fm)PM  =  AUkm  =  AmUkm  +  AwUkm  +  am, 

fk)Pyk)  -  fn)Py(m)  =  Avkm  =  Amvkln  +  A«vfa  +  A^y 


km 


(12) 


V(a)  =)P  + 


dw  ^  /  dz  jy 

v<r(fa)) 


(9) 


which  is  similar  to  McDougall’s  (1995)  formulation.  Here, 
yis  the  parameter  to  be  determined  by  the  thermal  wind 


Fig.  1.  Isoline  of  potential  vorticity  ( q )  is  the  trajectory  on  the 
isopycnal  surface. 
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Fig.  2.  Vertical  turning  of  the  P-vector:  (a)  P  vector  at  two 
isopycnal  levels,  and  (b)  turning  angle  between  two 
isopycnal  levels. 
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which  are  two  linear  algebraic  equations  for  yk)  and  fm). 
Here  y(l>  =  y(x,  y,  cr7),  and 


[A(0)ukm,A{0)vkm)  =  { 


)=r^*xvj 
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5.  P-Vector  Spiral 

McDougall  (1995)  first  presented  the  P-spiral  con¬ 
cept  on  the  neutral  surface.  Letting  <f>  be  the  angle  mea¬ 
sured  counterclockwise  from  the  x-axis  to  the  vector  P 
(McDougall,  1995),  we  have 


is  the  vertical  shear  of  y(<7)  with  no  ^-effect  and  no 
diapycnal  velocity.  With  the  /i-effect  and  diapycnal  ve¬ 
locity,  we  will  have  two  extra  source  terms  (A^ukm, 
A (wvAm)  and  (A (w)ukm,  AMvkm)  on  the  right-hand  sides  of 
(12).  As  soon  as  yk)  is  obtained,  the  velocity  field  can  be 
computed  by  (9)  or  (11). 

If  the  determinant 


cos  <j>  =  Px,  sin  0=P  (16) 

Existence  of  the  P-spiral  between  potential-density  lev¬ 
els  becomes  the  nonzero  d(f>/dz(<T) .  We  may  define  a 
nondimensional  parameter  (i.e.,  the  P-spiral  parameter) 
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the  algebraic  equations  (12)  have  definite  solutions  for 
yk)  (m  ±  k ): 
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The  determinant  in  the  denominator  of  (14)  is  the  sine  of 
the  vertical  turning  angle  between  Pik)  and  P‘m>  (Fig.  2(b)). 
The  solution  exists  if 


A  =  ^ 

3z(<7) 


(17) 


to  describe  the  P-spiral.  Substitution  of  (16)  into  (17)  leads 
to 


A  =  _/|(g)3Pv/3z(g)=/;(g)  dPy/dz 
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and  substitution  of  (10)  into  (18)  leads  to 
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where  <j>km  indicates  the  P-vector  turning  angle  between 
the  two  levels  ak  and  am. 

The  conventional  view  on  the  validity  of  the  inver¬ 
sion  is  the  existence  of  velocity  spiral,  which  is  called 
the  /1-spiral  by  Stommel  and  Schott  (1977).  This  is  true 
only  for  a  system  with  no  ^-effect  and  no  diapycnal  ve¬ 
locity.  For  a  system  with  either  the  ^-effect  or  the 
diapycnal  velocity,  the  validity  of  the  inversion  is  no 
longer  the  existence  of  velocity  spiral,  but  the  existence 
of  the  P-vector  spiral  ((/)km  ±  0).  Therefore,  we  have  the 
benefit  of  use  of  the  unit  vector,  P,  because  the  existence 
of  P  is  equivalent  to  the  intersection  of  isopycnals  and 
surfaces  of  constant  potential  vorticity  (necessary  condi¬ 
tion  1);  and  the  existence  of  the  P-spiral  is  required  for 
the  inversion  (necessary  condition  2). 


is  the  Jacobian. 

Thus,  the  necessary  condition  for  occurrence  of  the 
P-spiral  is  the  non-zero  Jacobian  of  qia)  and  dq,<y)/dz(c\ 
that  is,  the  nonzero  q(a)  and  d qia)ldz(a\  and  the  independ¬ 
ence  between  q[a)  and  3<5f(tTV3z( 

Substitution  of  (6)  into  (19)  leads  to 
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which  indicates  that  the  P-spiral  can  be  generated  by  lat¬ 
erally  and  vertically  inhomogeneous  thickness  between 
two  closely-spaced  potential-density  surfaces.  Recently, 
Chu  and  Li  (2000)  obtained  the  South  China  Sea  isopycnal 
surface  circulation  using  the  P-vector  sprial  method. 

6.  Climatological  Annual  Mean  Z(a>,  Q{a),  and  A  Data 

6.1  Division  of  potential  density  layers 

We  use  three  reference  levels  for  the  potential  den¬ 
sity  computation:  cr0  for  the  surface,  a2  for  2  km  depth, 
and  C74  for  4  km  depth.  Here, 

om  =  p-  1000  kg/m3,  m  =  0,2,4.  (21) 

The  potential  density  cr0,  cr2,  and  cr4  represent  upper,  in¬ 
termediate,  and  deep  layers,  respectively.  We  consider  the 
following  ranges  for  the  crffl-values: 

25.01  <  cr0<  27.99,  33.11  <  cr,  <  37.09, 

41.51  <ct4<  45.99.  (22) 

Notice  that  the  values  of  28.00,  37.10,  46.0  are  maximum 
values  for  c0,  c2,  and  cr4  computed  from  the  NODC  T,  S 
data  set.  We  further  discretize  c0,  cr2,  and  cr4  with  the  same 
increment, 

Act  =  0.02  kg/m3.  (23) 

Thus,  we  have  150  <70-layers,  200  oylayers,  and  225  c4- 
layers.  Within  each  layer,  the  density  is  vertically  uni¬ 
form. 


Pv.  ftk 


Fig.  3.  Vertical  discretization  of  the  multi-layer  ocean  with  the 
k-th  layer  having  potential  density  ph  layer  thickness  h  k, 
top  and  bottom  pressures  (  p  k,  Pkb),  respectively. 


refer  to  different  dimensions:  <7  (zf  for  the  depth  -z.j,  and 
a(k).  Thus,  we  use  the  sigma  value  at  the  bottom  of  the 
/r-th  isopycnal  layer,  cffk),  for  the  comparison, 

ah(k)  =  \[(j(k)  +  <j(k  +  l)].  (24) 

The  depth  for  the  bottom  of  the  a(k)- layer  is  obtained  by 

z[a)=Zj,  if  &{zj)  =  <Tb{k)  (25) 

and 


6.2  Establishment  of  hydrographic  data  with  high  verti¬ 
cal  resolution 

We  use  the  NODC  climatological  T,  S  data  for  the 
Atlantic  Ocean  in  this  study.  The  data  set  has  1°  x  1°  hori¬ 
zontal  resolution  and  33  vertical  levels  (Levitus  and 
Boyer,  1994;  Levitus  et  al.,  1994).  In  order  to  well  re¬ 
solve  potential-density  surfaces,  we  use  cubic  spline  to 
interpolate  the  T,  S  data  into  236  z-levels  with  three  dif¬ 
ferent  increments:  10  m  from  0  to  1000  m  depths,  20  m 
from  1000  to  2500  m  depths,  and  50  m  from  2500  to  5500 
m  depths.  Thus,  we  build  up  a  high  resolution  z-coordi- 
nate  data  set  [T{zf),  S(zj),  <70 (z,),  <72(z2),  <74(z;)].  The 
symbol  “A”  indicates  either  the  observational  data  or  the 
values  computed  from  the  observational  data. 

6.3  Transformation  of  T,  S  data  from  z-coordinates  to 
potential  density  surfaces 

The  transformation  is  fulfilled  by  comparing  the  z- 
coordinate  potential  density  data  <7  (zf  =  [  <7 0(z;),  ct2(zj), 
<74(z2)],  with  the  discrete  cr-values,  a(k).  The  two  values 


zia)  =  Zj  + 


o{k)~ 

-a( 

zi)  I 

<7| 

|-6 

iz){ 

(zDl  ~Zj\ 


if  a(zj)<<7h{k)<a(zj+1) 


(26) 


where  the  superscript  b  indicates  the  bottom  of  the  k-th 
isopycnal  layer.  The  thickness  of  the  k-th  isopycnal  layer 
is  obtained  by  (Fig.  3) 


For  the  k-th  potential-density  layer  <J(k),  the  poten¬ 
tial  vorticity  Q k(<7)  and  the  pseudo  potential  vorticity  qk(a) 
are  computed  by 


<r= = in 

hr. 


(27) 
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Fig.  4.  Climatological  (P ,  Q)  data  on  <J2  =  35.47  (kg/m3)  isopycnal  surface:  (a)  pressure,  P  (unit:  db),  and  (b)  potential  vorticity, 
Q  (unit:  10-n  m-1s-1). 


Having  the  climatological  ( hka\  Qka))  data  (Fig.  4),  it 
is  easy  to  invert  the  absolute  velocity  on  the  potential- 
density  surfaces.  After  checking  the  whole  data  set,  we 
found  the  existence  of  the  P-spiral  at  all  water  columns 
of  the  Atlantic  Ocean  on  1°  x  1°  grid.  Earlier  work  (e.g., 
McCartney,  1982;  Keffer,  1985;  Talley,  1988)  also  shows 
the  benefit  of  using  potential  vorticity  in  ocean  circula¬ 
tion  studies. 

6.4  Parameter  A 

The  spatial  variability  of  the  parameter  A  represents 
the  spatial  variability  of  the  P-vector  spirals.  For  illustra¬ 
tion,  we  chose  six  locations  (Fig.  5),  A  (54°W,  20°N),  B 
(36°W,  28°N),  C  (20°W,  55°N),  D  (57°W,  40°N),  E  (18°W, 
6°S),  and  F  (30°W,  45°N),  where  the  first  three  locations 
were  studied  extensively  in  the  past  (e.g.,  Schott  and 
Stommel,  1978;  Stommel  and  Schott,  1997).  These  six 
points  represent  the  Antilles  Current,  North  Atlantic  Gyre, 
Norwegian  Current,  Gulf  Stream,  South  Equatorial  Cur¬ 
rent,  Antarctic  Circumpolar  Current,  respectively.  We 
computed  the  parameter  A(ct2)  at  each  cr2-level  for  the 
six  locations  (Fig.  6).  The  values  of  A(ct2)  are  usually 
small  in  the  upper  isopycnal  layers,  and  become  large  in 
deeper  isopycnal  layers. 

To  quantitatively  estimate  the  turning  angle  between 
the  two  consecutive  isopycnal  layers,  we  discretize  (17) 
into 

A(f>  =  A.  (28) 

For  A  =  0.2,  the  vertical  turning  angle  is  around  11°. 
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Fig.  5.  Six  locations  chosen  for  illustration:  the  Antilles  Cur¬ 
rent  (A),  North  Atlantic  Gyre  (B),  Norwegian  Current  (C), 
Gulf  Stream  (D),  South  Equatorial  Current  (E),  Antarctic 
Circumpolar  Current  (F). 


7.  Zonally  Integrated  Latitudinal  Transport 

Following  Wunsch  and  Grant  (1982),  we  integrate 
the  inverted  velocity  v  zonally  (Fig.  7)  and  study  the  flow 
in  the  meridional  plane.  The  flow  (Fig.  7)  shows  similar 
patterns  as  depicted  in  Wunsch  and  Grant  (1982).  The 
dominant  features  are  poleward  movement  above  about 
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Fig.  6.  Climatological  P-vector  spiral  parameter  (A)  at  several  locations. 


(70  =  27.76  with  equatorward  movement  below  that,  in 
general  accord  with  the  conventional  picture  of  what 
should  happen.  The  maximum  northward  motion  in  the 
upper  layer  (cr0  >  27.76)  occurs  at  35°N  with  a  transport 
of  59.9  Sv  (1  Sv  =  106  m3/s).  Our  model  also  shows  the 
system  changes  from  a  net  downwelling  from  the  upper 
layers  into  the  deep  water  poleward  of  40°N,  to  a  net 
upwelling  through  the  water  column  equatorward  of  that 
latitude.  This  consists  with  earlier  studies  by  Wunsch  and 
Grant  (1982). 


8.  Conclusions 

(1)  In  this  study,  we  prove  that  the  conventional 
view  on  the  validity  of  the  inversion  (existence  of  veloc¬ 
ity  spiral)  is  true  only  for  a  system  with  no  /i-effect  and 
no  diapycnal  velocity.  For  a  system  with  either  the  /i-ef- 
fect  or  the  diapycnal  velocity,  the  validity  of  the  inver¬ 
sion  is  the  existence  of  the  P-spiral. 

(2)  The  P-spiral  occurs  when  the  potential  vorticity 
(q(a))  and  its  vertical  gradient  ( dqia)/dz(a) )  are  nonzero  and 
independent.  These  spirals  are  featured  by  the  rotation  of 
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Latitude  (N) 

Fig.  7.  Zonally  integrated  latitudinal  transport  (Sv)  between  two  adjacent  isopycnal  levels  inverted  from  NODC  climatological 
mean  T,  S  data.  The  size  of  the  vector  indicates  the  magnitude  of  the  volume  transport. 


the  unit  vector,  P ,  on  the  isopycnal  surface  with  potential 
density.  Using  the  climatological  T,  S  data,  we  found  that 
the  P-spirals  exist  almost  everywhere  in  the  world  oceans. 

(3)  Based  on  the  potential  vorticity  conservation, 
we  developed  the  P-spiral  method  to  invert  the  absolute 
geostrophic  velocity  on  the  potential-density  surfaces. 
This  method  can  invert  the  absolute  velocity  on  the  po¬ 
tential-density  surfaces  for  a  dynamical  system  includ¬ 
ing  the  /i-effect  and  the  diapycnal  velocity,  and  has  great 
potential  for  investigating  the  general  ocean  circulations. 
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